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Abstract 


This  report  presents  models  of  seaways  that  can  be  used  for  simulation  of  ship  mo¬ 
tions.  The  present  work  assumes  deep  water,  in  which  ocean  wavelengths  are  less 
than  half  the  water  depth.  Seaway  models  are  presented  for  both  regular  and  ran¬ 
dom  wave  conditions,  and  are  given  in  both  fixed  and  translating  axis  systems.  For 
regular  waves,  the  nonlinear  free  surface  can  be  modelled  using  Stokes  second  order 
theory.  For  both  regular  and  random  waves,  the  influence  of  the  location  of  the  free 
surface  on  wave  kinematics  can  be  modelled  using  Wheeler  stretching.  Commonly 
used  wave  spectra  and  sources  of  wave  climate  data  are  presented  as  reference  in¬ 
formation.  Sample  computations  for  the  frigate  HALIFAX  illustrate  the  influence 
of  simulation  parameters  on  predicted  RMS  motions  and  required  computational 
time.  Required  computational  time  increases  significantly  when  including  nonlin¬ 
ear  forces  from  buoyancy  and  incident  waves.  In  most  cases,  simulations  run  faster 
than  real  time. 


Resume 


Dans  ce  rapport,  nous  presentons  des  modeles  de  houle  pouvant  servir  a  simu- 
ler  les  mouvements  de  navires.  Dans  cet  article,  nous  supposons  des  conditions 
d’eaux  profondes  dans  lesquelles  la  longueur  d’onde  des  vagues  est  inferieures  a  la 
moitie  de  la  profondeur.  Nous  presentons  des  modeles  de  houle  pour  des  condi¬ 
tions  de  vagues  regulieres  et  aleatoires,  et  ces  modeles  sont  donnees  pour  des 
systemes  de  coordonnees  fixes  et  en  translation.  Dans  le  cas  de  vagues  regulieres, 
on  peut  modeliser  la  surface  libre  non  lineaire  en  appliquant  la  theorie  de  Stokes  de 
deuxieme  ordre.  Dans  le  cas  de  vagues  regulieres  et  aleatoires,  on  peut  modeliser 
l’effet  de  l’emplacement  de  la  surface  libre  sur  la  cinematique  des  vagues  en  utili- 
sant  un  etalement  de  Wheeler.  Nous  presentons  a  titre  d’information  de  reference 
les  spectres  de  vagues  couramment  utilises  et  les  sources  de  donnees  sur  les  vagues. 
Des  calculs  echantillons  pour  la  fregate  HALIFAX  illustrent  l’effet  des  parametres 
des  simulations  sur  la  valeur  quadratique  des  mouvements  prevus  et  sur  le  temps 
de  calcul  necessaire.  Le  temps  de  calcul  necessaire  augmente  considerablement 
lorsqu’on  inclut  les  forces  non  lineaires  associees  a  la  flottabilite  et  aux  vagues  in- 
cidentes.  Dans  la  plupart  des  cas,  les  simulations  sont  plus  rapides  que  le  processus 
simule  en  temps  reel. 
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Executive  summary 

Introduction 

DRDC  Atlantic  is  developing  a  new  object-oriented  library  for  simulation  of  ship 
motions  in  waves.  This  report  presents  models  of  regular  and  random  seaways  that 
can  be  used  during  simulation  of  ship  motions.  The  present  work  assumes  deep 
water,  in  which  seaway  wavelengths  are  less  than  half  the  water  depth. 

Principal  Results 

Models  of  regular  and  random  seaways  have  been  developed  using  both  fixed  and 
translating  coordinate  systems.  The  relationship  between  wave  phases  for  fixed  and 
translating  coordinate  systems  is  presented.  For  regular  waves,  the  nonlinear  free 
surface  can  be  modelled  using  Stokes  second  order  theory.  For  both  regular  and 
random  waves,  the  influence  of  the  location  of  the  free  surface  on  wave  kinematics 
can  be  modelled  using  Wheeler  stretching.  Sample  motion  computations  for  the 
frigate  HALIFAX  indicate  that  required  computational  time  increases  significantly 
when  including  nonlinear  forces  from  buoyancy  and  incident  waves.  In  most  cases, 
inclusion  of  nonlinear  buoyancy  and  incident  wave  forces  has  little  influence  on 
predicted  RMS  motions. 

Significance  of  Results 

Seaway  models  are  now  available  for  simulation  applications  such  as  predicting 
ship  motions  and  visualizing  seaways.  The  presented  wave  phase  relationship  be¬ 
tween  fixed  and  translating  axis  systems  can  be  used  for  applications  such  as  pre¬ 
dicting  motions  of  a  ship  freely  maneuvering  in  a  seaway  or  predicting  motions  of 
several  ships  in  a  seaway.  Presented  ship  motion  predictions  typically  run  faster 
than  real  time,  indicating  they  are  suitable  for  real  time  simulation  applications. 


Future  Plans 

Development  of  DRDC  Atlantic’s  object-oriented  ship  motion  library  is  continuing. 
The  next  phase  of  work  will  incorporate  forces  from  rudders,  bilge  keels,  and  other 
appendages  into  ship  motion  predictions. 
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Sommaire 


Introduction 

RDDC  Atlantique  est  en  train  d’elaborer  une  nouvelle  base  de  donnees  orientee  objet  pour 
la  simulation  des  mouvements  des  navires  dans  la  houle.  Dans  ce  rapport,  on  presente  des 
modeles  de  houle  reguliere  et  aleatoire  pouvant  servir  lors  de  la  simulation  des  mouve¬ 
ments  des  navires.  Dans  la  presente  etude,  on  suppose  des  eaux  profondes  dans  lesquelles 
la  longueur  d’onde  des  vagues  est  inferieure  a  la  moitie  de  la  profondeur. 

Resultats  principaux 

Nous  avons  elabore  des  modeles  de  houle  reguliere  et  aleatoire,  en  utilisant  des  systemes  de 
coordonnees  fixes  et  en  translation.  Nous  presentons  la  relation  entre  les  phases  des  vagues 
dans  ces  deux  systemes  de  coordonnees.  Des  calculs  echantillons  de  mouvements  pour  la 
fregate  HALIFAX  indiquent  que  le  temps  de  calcul  necessaire  augmente  considerablement 
lorsqu’on  inclut  les  forces  non  lineaires  associees  a  la  flottabilite  et  aux  vagues  incidentes. 
Dans  la  plupart  des  cas,  l'inclusion  de  ces  forces  non  lineaires  n’a  guere  d’effet  sur  la  valeur 
quadratique  des  mouvements. 

Importance  des  resultats 

Nous  disposons  maintenant  de  modeles  de  houle  qui  permettent  de  proceder  a  des  simula¬ 
tions  comme  la  prevision  des  mouvements  d’un  navire  et  la  visualisation  de  la  houle.  La 
relation,  presentee  dans  cet  article,  entre  les  systemes  de  coordonnees  fixes  et  en  translation 
qui  decrit  les  phases  des  vagues  peut  etre  utilisee,  par  exemple,  pour  prevoir  les  mouvement 
de  navires  qui  manceuvrent  librement  dans  la  houle  ou  encore  pour  prevoir  les  mouvements 
de  plusieurs  navires  dans  la  houle.  Les  previsions  de  mouvement,  que  nous  presentons 
dans  cet  article,  precedent  normalement  plus  rapidement  qu’en  temps  reel,  ce  qui  indique 
qu’elles  conviennent  aux  simulations  en  temps  reel. 

Travaux  ulterieurs  prevus 

L’ elaboration  d’une  nouvelle  base  de  donnees  orientee  objet  sur  les  mouvements  des  navires 
dans  la  houle  se  poursuit.  Dans  la  prochaine  phase  des  travaux,  nous  inclurons  dans  les 
previsions  les  forces  exercees  sur  le  gouvernail,  la  carene  et  les  autres  surfaces. 
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1  Introduction 


DRDC  Atlantic  is  developing  software  components  for  modelling  and  simulation  of 
ships  in  waves.  Work  completed  to  date  includes  predictions  of  forces  and  motions 
in  waves  for  an  unappended  ship  hull  in  the  frequency  domain  [1]  and  time  domain 
[2]. 

This  report  presents  models  of  seaways  that  can  be  used  for  various  simulation  ap¬ 
plications,  including  prediction  of  ship  motions  and  sea  loads.  Models  of  regular 
seaways  are  useful  for  comparison  of  numerical  predictions  with  experimental  re¬ 
sults,  and  for  computations  in  the  frequency  domain.  Models  of  random  seaways 
are  useful  for  modelling  realistic  sea  conditions. 

When  computing  ship  motions  in  waves,  it  is  often  convenient  to  use  a  translating 
coordinate  system  that  moves  with  the  forward  speed  of  the  ship.  This  report  con¬ 
siders  seaway  modelling  in  both  fixed  and  translating  coordinate  systems,  including 
wave  phasing  relationships  between  the  two.  Correct  treatment  of  wave  phasing  is 
necessary  for  simulation  applications  such  as  hull  slamming  and  motions  of  two 
ships  during  replenishment  at  sea. 

The  present  work  is  limited  to  deep  seaways,  in  which  the  presence  of  the  ocean 
bottom  has  no  influence  on  water  waves.  In  practical  terms,  a  seaway  can  be  con¬ 
sidered  deep  when  the  water  depth  is  greater  than  half  the  length  of  incident  waves. 
Most  ship  motion  predictions,  including  those  in  References  1  and  2,  are  based  on 
the  assumption  of  deep  seaways. 

Section  2  describes  the  treatment  of  regular  seaways  in  both  the  frequency  and 
time  domains.  Simulation  of  realistic  random  seaways  is  described  in  Section  3, 
followed  by  descriptions  of  commonly  used  wave  spectral  models  in  Section  4. 
Section  5  presents  background  information  on  wave  climate  data  relevant  for  devel¬ 
oping  simulations.  Section  6  presents  an  improved  treatment  of  rotational  moments 
arising  from  nonlinear  forces  due  to  buoyancy  and  incident  waves.  Implementation 
of  seaway  models  using  the  Python  programming  language  is  described  in  Sec¬ 
tion  7.  Section  8  describes  computational  requirements  for  simulating  motions  of  a 
frigate  in  a  seaway,  and  is  followed  by  comparisons  between  experimental  and  pre¬ 
dicted  motions  of  a  frigate  in  Section  9.  Final  conclusions  are  given  in  Section  10. 

2  Regular  Seaways 


Seaways  comprised  of  regular  waves  (i.e.,  waves  of  a  single  frequency  and  direc¬ 
tion)  are  of  fundamental  importance  for  ship  motion  predictions.  Motions  in  ran¬ 
dom  seaways  can  be  evaluated  using  linear  superposition  of  motion  predictions 
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from  regular  waves.  Due  to  their  repeatability,  regular  seaways  are  also  important 
for  validation  of  numerical  ship  motion  predictions  using  experiments. 

Several  texts,  including  Sarpkaya  and  Isaacson  [3]  and  Chakrabarti  [4],  give  de¬ 
tailed  information  regarding  the  theory  of  regular  waves.  The  present  report  as¬ 
sumes  that  the  water  is  sufficiently  deep  that  the  ocean  bottom  does  not  influence 
the  waves.  In  practice,  bottom  effects  are  negligible  when  the  water  depth  is  greater 
than  half  the  wavelength. 

Regular  waves  can  be  considered  in  various  different  coordinate  systems.  It  is  sim¬ 
plest  to  consider  waves  in  a  fixed  coordinate  system.  A  fixed  coordinate  system 
is  suitable  for  evaluating  loads  on  fixed  structures  such  as  stationary  offshore  plat¬ 
forms.  When  considering  wave  forces  acting  on  a  moving  object,  such  as  a  ship 
with  forward  speed,  it  is  more  appropriate  to  treat  waves  using  a  moving  coordi¬ 
nate  system.  The  present  section  begins  with  a  treatment  using  fixed  coordinates 
followed  by  the  more  complicated  case  of  a  translating  coordinate  system. 


2.1  Linear  Regular  Seaway  in  Fixed  Coordinates 

A  regular  seaway  is  most  frequently  modelled  using  sinusoidal  waves.  The  appli¬ 
cation  of  sinusoidal  waves  to  a  seaway  is  based  on  the  assumption  that  the  wave 
amplitude  is  small  relative  to  the  wavelength.  In  the  present  discussion,  incident 
waves  alone  are  considered,  defined  as  those  waves  which  are  not  influenced  by  the 
presence  of  a  body  in  the  seaway.  In  reality,  the  presence  of  a  body  such  as  a  ship 
or  offshore  structure  in  a  seaway  will  produce  diffracted  waves  in  the  vicinity  of  the 
body. 


Figure  1  shows  the  fixed  coordinate  system  used  for  evaluating  waves.  The  waves 
are  travelling  from  direction  u.  The  incident  wave  elevation  in  fixed  coordinates  is 
given  by: 


C i(xf,yf,t )  = 


a  cos 


hr  ( y] 


sm  u 


x *  cos  v 


)  -ix it 


(1) 


where  a  is  the  wave  amplitude,  kj  is  the  incident  wavenumber,  x j  and  yj  are  the 
horizontal  plane  coordinates,  c oi  is  the  incident  wave  frequency,  and  is  the  phase 
lead  angle  for  the  wave  crest  at  the  origin.  In  the  present  work  assuming  deep  water, 
the  following  dispersion  relation  applies: 

hT  =  —  (2) 

9 

where  g  is  gravitational  acceleration.  The  velocity  potential  and  its  derivatives  are: 


_  9  a  ■ 

<h/  =  —  sm 

Ui 


ki  (y*  sin  v  —  x*  cos  v)  —  ui  t 


exp  (kj  zw{)  (3) 
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Ship 


Figure  1:  Fixed  Coordinate  System 


<9$/ 

dxf 

=  —a  ujj  cos  v  cos  kj  (y*  sin  v  —  cos  u)  —  loj  t 

-  ef] 

ti 

x  exp (&/  zwi) 

(4) 

dyf 

=  a  uji  sin  v  cos  kj  {]/  sin  v  —  x*  cos  v)  —  u>i  t  - 

-  ef] 

ti 

x  exp (&/  zwl) 

(5) 

dz 

=  a  LUi  sin  kj  ( y 7  sin  v  —  x^  cos  v)  —  i ojt  — 

exp  [kj  zwi ) 

(6) 

dt 

=  —g  a  cos  kj  (yf  sin  v  —  x?  cos  v)  —  Uj  t  — 

exp  (fc/  zw{) 

(7) 

where  zwi  is  the  elevation  (positive  upward)  relative  to  the  calm  water  surface. 
When  applying  the  above  equations,  it  should  be  noted  that  the  velocity  potential 
and  its  derivatives  will  be  zero  for  locations  above  the  free  water  surface. 

2.2  Second  Order  Stokes  Theory  for  Regular 
Seaway  in  Fixed  Coordinates 

Various  nonlinear  models  are  available  which  provide  improved  accuracy  for  regu¬ 
lar  seaways,  as  discussed  in  References  3,  4  and  5.  The  most  widely  used  nonlinear 
wave  theory  is  Stake’s  second  order  theory.  In  deep  water,  the  wave  elevation  given 
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Stokes  second  order  theory 
Linear  theory 


Figure  2:  Profiles  of  Regular  Wave  at  Limiting  Steepness  H / A  =  0.142  (no 
distortion  of  axes) 


by  Stoke ’s  second  order  theory  is: 


(i(x*  ,y-f  ,t)  =  a  cos  kj  (y*  sinz/  —  x*  cost/)  —  cu/ f  — 


a 2  k 


cos  2 


ki  (y1 


sin  u 


_  q-f 


XJ  COS  V 


e{ 


(8) 


For  second  order  Stokes  waves  in  deep  water,  the  velocity  potential  and  its  deriva¬ 
tives  are  the  same  as  for  linear  theory  given  in  Equations  3  to  7.  Similarly  the 
dispersion  relation  for  Stokes  second  order  theory  is  the  same  as  that  for  linear 
theory  given  in  Equation  (2). 


Due  to  wave  breaking,  the  maximum  possible  wave  steepness  H / A  in  deep  water 
is  0.142,  where  H  is  wave  height  and  A  is  wavelength.  Figure  2  gives  wave  profiles 
from  linear  and  second  order  theories  at  this  limiting  wave  steepness.  At  lower 
wave  steepnesses,  the  difference  between  linear  and  second  order  profiles  will  be 
less  pronounced. 


2.3  Regular  Seaway  in  Translating  Earth 
Coordinates 

For  a  ship  travelling  through  a  seaway,  the  forward  speed  of  the  ship  has  a  signifi¬ 
cant  influence  on  the  interactions  between  the  ship  and  the  seaway.  For  example,  the 
frequency  of  waves  encountered  by  the  ship  will  differ  from  the  incident  wave  fre¬ 
quency  if  the  ship  is  travelling  with  forward  speed  at  any  heading  other  than  beam 
seas.  To  account  for  the  influence  of  ship  forward  speed,  ship  motion  computations 
are  routinely  performed  in  translating  earth  coordinates.  The  ship  is  assumed  to 
travel  at  steady  forward  speed  and  heading,  and  the  coordinate  system  translates 
accordingly.  Figure  3  shows  the  translating  earth  coordinate  system  and  Figure  4 
shows  the  sea  direction  relative  to  ship  forward  speed.  The  origin  of  the  translating 
earth  coordinate  system  is  located  at  the  ship  centre  of  gravity  for  the  ship  in  calm 
water.  When  evaluating  wave  kinematics,  the  vertical  coordinate  zwi  relative  to  the 
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calm  waterline  is  used.  The  following  equation  relates  the  relative  sea  direction  to 
terms  of  the  fixed  axis  coordinate  system  presented  in  Figure  1 . 

ft  =  u  +  180°  -  X  (9) 

where  x  is  the  ship  heading  relative  to  the  x?  axis. 

For  linear  regular  waves,  the  wave  elevation  in  translating  earth  axes  is  given  by: 

C i(x,y,t )  =  a  cos  [kj  (x  cos  ft  —  y  sin/3)  —  c oet  —  eft  (10) 

where  x  and  y  are  horizontal  plane  coordinates  as  shown  in  Figure  3,  ue  is  the  wave 
encounter  frequency,  and  ej  is  the  wave  phase  lead  of  the  crest  at  the  origin.  The 
wave  encounter  frequency  is  given  by: 

c oe  =  \uji  —  U  kj  cos  ft\  (11) 

where  U  is  the  ship  forward  speed.  The  wave  velocity  potential  and  its  derivatives 
are  given  by: 


<9$j 

dx 

<9$7 

dy 

0$/ 

dz 

~dt 


- —  sin  [kj  (x  cos/3  —  y  sin/3)  —  c oet  —  eft  exp  (ki  zwft)  (12) 

Ui 

a  Uj  cos  ft  cos  [kj  (x  cos  ft  —  y  sin  ft)  —  u>et  —  eft 
x  exp (fc7  zwi)  (13) 

— a  ui  sin/3  cos  [kj  (xcos/3  —  y  sin  ft)  —  uet  —  eft 
x  exp (fc7  zwi)  (14) 

a  ui  sin  [kj  {x  cos  ft  —  y  sin/3)  —  c oet  —  eft  exp  [ki  zwft  (15) 


—  - — uje  cos  [ki  (x  cos  ft  —  y  sin/3)  —  c oet  —  eft 

ui 

x  exp (kT  zwi) 


(16) 


If  the  regular  seaway  is  modelled  using  Stokes  second  order  theory,  then  the  wave 
elevation  in  translating  earth  axes  is: 

(i(x,y,t)  =  a  cos  [kj  (x  cos  ft  —  y  sin/3)  —  uet  —  eft 
a2ki 

4 - —  cos  2  [hi  (x  cos  ft  —  y  sin/3)  —  uet  —  eft  (17) 


2.4  Frequency  Domain  Representation  of  Regular 
Seaway  in  Translating  Earth  Coordinates 


Ship  motions  are  often  computed  in  the  frequency  domain  using  translating  earth 
coordinates,  as  described  in  References  1  and  2.  Relationships  between  terms  in 
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the  time  and  frequency  domains  can  be  expressed  by: 


(18) 

(19) 

(20) 


r)j(t)  =  Real  {r)j(ue)  exp(i  toe  t)}  for  j  =  1  —  6 
C (t)  =  Real{C(cce)  exp(i  tve  t)j 
<f>/(f)  =  Real{0/(u;e)  exp(i  cze  t)} 

where  rj3  is  ship  motion  displacement  for  mode  j.  The  frequency  domain  terms 
r)j( u>e),  ((u-'p),  and  <j>i(ue)  in  the  above  equations  are  complex  quantities  which  in¬ 
clude  both  amplitude  and  phase  information.  For  the  incident  wave  velocity  poten¬ 
tials  given  in  Equation  (20),  further  distinction  between  time  domain  and  frequency 
domain  quantities  is  made  by  using  upper  and  lower  case  letters  $/(£)  and  0/(u/e). 

Incident  wave  elevation  in  the  frequency  domain  is  given  by: 

Ci  —  a  exp  [—ikj  (a;  cos  (3  —  y  sin  /3)  +  i  e{\  (21) 

The  incident  wave  potential  and  its  derivatives  are  given  in  the  frequency  domain 
by: 

— —  exp  [— %  ki  (x  cos (3  —  y  sin/3)  +  i  ej]  exp (ki  zwi)  (22) 

Ui  ~ 

a  u>i  cos  (3  exp  [—ikj  ( xcos/3  —  ysin/3)  +  i  ej\  exp  (kj  zwi)  (23) 
— a  ui  sin/3  exp  [— i  kj  (xcos/3  —  ysin/3 )  +  it/]  exp  (ki  zwi)  (24) 
i  a  uj  exp  [— %  kj  (x  cos  (3  —  y  sin/3)  +  i  ej\  exp  (ki  zwi)  (25) 

- — —  exp  [— i  kj  (x  cos  (3  —  //sin/3)  +  i  e  A  exp  (ki  zwA  (26) 

ui 

2.5  Relationships  Between  Fixed  and  Translating 
Earth  Coordinate  Systems 

Within  modelling  and  simulation  applications,  it  is  often  useful  to  be  able  to  shift 
from  a  fixed  to  a  translating  earth  coordinate  system.  One  example  is  a  ship  maneu¬ 
vering  in  a  seaway,  whereby  many  wave-induced  force  terms  are  most  conveniently 
evaluated  by  transformation  to  translating  earth  coordinate  systems  at  each  time 
step.  If  the  origin  of  the  translating  earth  coordinate  system  is  located  at  x?,  yJ  at 
time  to,  then  the  phase  lead  in  translating  earth  axes  will  be  given  by: 

e/  =  +  to  (u)i  —  cue)  —  kj  (yf  sin  v  —  cos  v)  (27) 


07  = 

<90/ 

dx 
<907 
<9  y 

<907  _ 

dz 
<907 
<9 1 
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If  the  ship  is  assumed  to  maintain  steady  ship  speed  U  and  heading  u,  then  the 
following  relationships  will  hold: 


xf  (t) 

=  xf(t0)  +  U  (t 

—  to)  COS  V 

(28) 

yf(t ) 

=  yf(t0 )  -  U  {t 

—  to)  sin  v 

(29) 

<9$ 

d<l> 

(30) 

dx 

=  o  f  cos  x 

OX 1 

sm  x 

dy* 

<9$ 

d<l> 

(31) 

dy 

=  cos  x  + 

oyf 

sm  x 

ox J 

In  the  context  of  computing  wave  forces  during  simulation  of  a  ship  in  seaway, 
the  assumption  of  steady  speed  and  heading  will  have  practical  implications.  If 
the  ship  is  assumed  to  have  steady  speed  and  heading,  the  seaway  will  require 
transformation  of  the  seaway  from  earth-fixed  axes  to  translating  earth  axes  only 
at  the  beginning  of  the  simulation.  Alternatively,  if  a  simulation  models  a  freely 
maneuvering  ship  with  speed  and  heading  varying  at  each  time  step,  then  wave 
force  calculations  will  require  transformation  of  the  seaway  from  earth-fixed  axes 
to  translating  earth  axes  at  each  time  step. 


2.6  Modelling  of  Finite  Amplitude  Wave  Kinematics 
Using  Wheeler  Stretching 

The  wave  kinematic  terms,  presented  in  Equations  3  to  7  for  fixed  axes  and  in  Equa¬ 
tions  12  to  16  for  translating  axes,  are  based  on  the  assumption  of  small  amplitude 
waves.  As  wave  amplitude  becomes  increasingly  large,  questions  arise  regarding 
the  validity  of  the  equations  for  wave  kinematic  terms,  particularly  in  the  vicinity 
of  the  instantaneous  water  surface.  Xu  [6]  compares  various  models  used  for  wave 
kinematic  terms  in  finite  amplitude  waves.  Among  the  models  discussed,  Wheeler 
stretching  [7]  has  been  selected  here  because  of  its  relative  simplicity  and  good 
agreement  with  more  complicated  wave  models.  In  Wheeler  stretching,  the  term 
exp  kizwi  is  replaced  by  exp  ki  (zwi  —  £(£)),  with  ((t)  being  the  instantaneous 
waterline. 

3  Simulation  of  Realistic  Random  Seaways 


Discussion  thus  far  has  been  limited  to  regular  waves,  limited  to  a  single  wave 
frequency  and  direction.  Realistic  seaways  consist  of  wave  components  arising 
from  many  frequencies  and  directions. 
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3.1  Random  Unidirectional  Seaway 

Random  seaways  can  be  simulated  by  linear  superposition  of  a  finite  number  of  reg¬ 
ular  wave  components.  Chakrabarti  [4]  provides  a  useful  overview  of  this  approach. 
Using  the  fixed  axis  systems  of  Figure  1,  the  wave  elevation  for  a  unidirectional 
random  seaway  can  be  simulated  as  follows: 

JVj 

c i(xf,yf,t)  =  Yai  cos 

i=  1 

where  Nj  is  the  number  of  regular  wave  components  used  to  represent  the  random 
seaway,  a*  is  the  wave  amplitude  for  component  i,  ki-i  is  the  wavenumber  for 
component  i,  ujj_t  is  the  incident  wave  frequency  for  component  i,  and  _t  is  the 
phase  for  wave  component  i.  The  wave  amplitude  for  component  i  is  based  on  a 
specified  wave  spectrum  as  follows: 

di  =  A(u u-i)  (33) 

where  SUl  is  the  spectral  density  and  A(ujr_t)  is  the  wave  frequency  interval  for 
component  i.  The  wave  phase  components  ej_{  are  obtained  using  random  number 
generation  based  on  a  uniform  distribution  between  0  and  2tt.  The  following  rela¬ 
tionship  exists  between  the  wave  component  amplitudes  and  the  standard  deviation 
of  water  elevation: 

1  Nl 

a2(0  =  2  £a*2  (34) 
1=1 

The  fidelity  of  the  simulated  seaway  will  increase  with  number  of  wave  compo¬ 
nents.  A  minimum  of  20  components  is  typically  used  for  simulating  unidirectional 
random  seaways. 

Wave  kinematic  terms  for  a  random  seaway  are  easily  determined  using  summation 
from  regular  wave  components.  When  applying  Wheeler  stretching  to  wave  kine¬ 
matic  terms  (Section  2.6),  the  instantaneous  elevation  should  be  evaluated  using 
contributions  from  all  wave  components. 


ki-i  (V 


sm  v 


x *  cos  v 


)  & I—i  t 


- I-i 


(32) 


3.2  Random  Directional  Seaway 


The  simulation  of  a  random  directional  seaway  includes  wave  components  from 
multiple  directions.  The  following  equation  is  used  for  modelling  the  wave  surface 
for  a  random  directional  seaway: 


JVj 


C  i(xf,yf,t)  = 


CLi  cos 


i= 0 


ki-i  ( y 1 


sin  Vi  —  x*  cos  Vi 


i)  -  Ui-i  t  - 


(35) 
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where  is  the  direction  for  wave  component  i.  The  amplitude  of  each  wave  com¬ 
ponent  is  given  by: 

dj  \J 2  f  u  i/jj j —j..  Uj )  A(cu/^i)  A  (r'j)  (36) 

where  SUltU  (cc/_j,  z/; )  is  the  directional  spectral  density  and  A(z/j)  is  the  direction 
interval  associated  with  wave  component  i. 

3.3  Nonlinear  Effects  for  Random  Seaways 

When  modelling  random  seaways,  the  only  nonlinear  effect  currently  modelled  is 
Wheeler  stretching,  which  is  applied  in  the  same  manner  as  for  regular  seas  (Sec¬ 
tion  2.6).  Forristall  [8]  gives  an  excellent  overview  of  methods  for  including  non¬ 
linear  effects  in  random  seaways.  It  appears  that  more  numerically  intensive  treat¬ 
ments  do  not  give  significantly  better  results  than  Wheeler  stretching. 

4  Commonly  Used  Wave  Spectral  Models 

The  simulation  of  random  seaways  described  in  the  previous  section  requires  input 
wave  spectral  densities.  Wave  spectra  are  typically  represented  using  commonly 
available  models,  such  as  those  presented  in  the  SHIPM07  manual  [9].  This  section 
gives  some  commonly  used  wave  spectral  models. 

4.1  Point  Wave  Spectra 

Point  wave  spectra,  denoted  SUl  (cu/),  give  wave  spectral  energy  as  a  function  of 
wave  frequency,  and  do  not  include  directional  information.  Point  wave  spectra  are 
often  applied  with  the  assumption  that  the  seaway  is  unidirectional  (i.e.,  all  waves 
are  travelling  from  the  same  direction).  This  assumption  is  typically  more  valid  in 
higher  sea  states,  while  lower  sea  states  are  typically  more  directionally  confused. 

The  RMS  (root-mean-square)  wave  elevation  is  related  to  the  wave  spectrum  by: 

o-(C)  =  y/rno  (37) 

/»oo 

m.i  =  /  ^(cu/)  dec/  (38) 

Jo 

where  m*  is  the  i’th  moment  for  the  wave  spectrum.  A  seaway  is  commonly  de¬ 
scribed  in  terms  of  significant  wave  height  Hs,  which  is  the  average  height  of  the 
highest  one  third  of  the  waves.  Assuming  a  Rayleigh  distribution  for  wave  eleva¬ 
tion,  significant  wave  height  is  given  by: 

Hs  =  4.005  y/m^  (39) 
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The  frequency  content  of  a  seaway  is  commonly  described  using  a  characteris¬ 
tic  wave  period,  for  which  several  different  definitions  are  available.  The  peak 
wave  period  Tp  is  the  wave  period  associated  with  the  wave  frequency  c uj  at  which 
Su,j(u)i)  has  its  maximum.  The  average  wave  period  Ti  is  another  characteristic 
wave  period,  and  is  given  by: 


T\ 


2  7T 


m  o 

m  i 


(40) 


A  third  commonly  used  period  is  the  zero-crossing  period  given  by: 


Tx 


(41) 


4.1.1  Bretschneider  Spectrum 


The  Bretschneider  spectrum  is  the  most  commonly  used  model  of  wave  spectra  in 
the  open  ocean.  Based  on  the  15th  International  Towing  Tank  Conference  (ITTC) 
[10],  the  formulation  for  the  Bretschneider  spectrum  is: 


SWI  (ui) 


486.0  Hg 


exp 


—  1948.2’ 


(42) 


The  above  spectrum  is  defined  in  terms  of  peak  wave  period  Tp.  For  a  Bretschneider 
spectrum,  the  following  relations  exist  with  the  average  and  zero-crossing  wave 
periods: 


Ti  =  0.773  Tp  (43) 

Tz  =  0.710  Tp  (44) 


4.1 .2  Three  Parameter  JONSWAP  Spectrum 


The  JONSWAP  spectrum  models  relatively  high-peaked  spectra  typically  encoun¬ 
tered  in  fetch-limited  regions  [10].  The  JONSWAP  spectrum  is  obtained  by  mul¬ 
tiplying  the  Bretschneider  spectrum  by  a  peak  enhancement  factor  accounting  for 
fetch-limited  conditions,  giving  the  following  [4]: 


SUl{ui) 


hi 


a 


a 


,  ,4 

jj2  ^ p 

44  s  5 


exp 


LUp 

-1.25— 7 

OJ4 


7 


exp 


(cv  -  Up)2 
2<72lu% 

fora;  <  up 
for  lu  >  up 


(45) 

(46) 

(47) 
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where  up  is  the  peak  wave  frequency  and  7  is  an  input  spectral  peak  parameter. 
Goda  [11]  derived  the  following  approximate  expression  for  the  normalization  term 


a 


0.0624 


a  = 


0.230  +  0.0336  7  -  0.185/(1.9  +  7) 


(48) 


The  JONS  WAP  spectrum  is  often  presented  as  a  two  parameter  spectrum,  with  the 
spectral  peak  parameter  7  having  a  default  value  of  3.3. 


4.1.3  Ochi  and  Hubble  Six  Parameter  Spectrum 


The  Ochi  and  Hubble  6  parameter  spectrum  [12]  models  collinear  swell  and  sea 
components  as  follows: 


UJl  (+f )  —  ^ 


2 

£ 

i= 1 


^P-i]Xx  h2s-i  exP  -  (- £ (^) 


r(Aj) 


(49) 


where  A*,  /is_i5  and  ccp_j  are  the  spectra  shape  parameter,  significant  wave  height, 
and  peak  frequency  for  component  i.  The  term  TfA/)  is  the  Gamma  function  with 
argument  A,  .  If  only  one  of  the  two  components  is  considered  and  the  shape  param¬ 
eter  A i  equals  one,  then  the  six  parameter  spectrum  is  equivalent  to  the  Bretschnei- 
der  spectrum. 


4.2  Directional  Wave  Spectra 

The  point  wave  spectra  discussed  above  do  not  include  directional  information.  To 
model  directionally  confused  seas,  directional  wave  spectra  denoted  SUuV  (u 7,  v) 
can  be  used. 


4.2.1  Point  Spectrum  with  Cosine  Squared  Spreading  Function 

A  directional  wave  spectrum  can  be  most  easily  modelled  by  multiplying  a  point 
spectrum  by  a  directional  spreading  function  as  follows: 

SUltV(u}i,v)  =  Sujiu^Giv)  (50) 

where  G{y)  is  a  directional  spreading  function.  Note  that  the  area  under  the  di¬ 
rectional  spreading  function  must  be  unity.  A  cosine  squared  spreading  function  is 
commonly  used,  with  the  following  format: 


G(u) 

1  2 

=  —  cos 

Is 

/  V  —  V  7rX 

V  7s  2J 

for  \u  —  v\  <  7S 

(51) 

G{6) 

=  0 

for  z/  —  v\  >  7 s 

(52) 
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where  v  is  the  principal  wave  direction  and  7.,  is  the  spreading  angle  in  radians.  A 
spreading  angle  of  7t/2  (90  degrees)  is  often  used  for  seakeeping  computations. 


4.2.2  Ten  Parameter  Directional  Spectrum 

Directional  seas  can  be  most  apparent  when  sea  and  swell  components  are  simi¬ 
lar  in  magnitude  and  are  approaching  from  different  directions.  The  ten  parameter 
spectrum  developed  by  Hogben  and  Cobb  [13]  is  suitable  for  modelling  such  con¬ 
ditions,  and  has  been  applied  to  ship  motion  predictions  by  Graham  and  Juszko 
[14,  15].  The  ten  parameter  spectrum  is  a  directional  extension  of  the  Ochi  and 
Hubble  six  parameter  spectrum,  with  each  of  the  swell  and  sea  components  being 
multiplied  by  its  own  directional  spreading  function  as  follows: 


Mi(v)  =  A(Pi)  cos2Pi  y  2  U,j  for  i  —  1, 2  and  \u  —  1 7)  <  vr/2  (53) 

Miiy)  =  0  for  %  =  1, 2  and  \v  —  z/j|  >  7t/2  (54) 

where  Pt  and  1 7  are  the  directional  spread  parameter  and  mean  direction  (from)  for 
component  i.  The  normalization  factor  A ( P, )  is  expressed  as: 


m) 


2(2P,-1)  p2 (p.  + 
7T  T(2Pi  +  1) 


i  =  1,2 


(55) 


where  T(2 Pi  +  1)  is  the  Gamma  function  with  argument  2 1\  +  1. 


5  Wave  Climate  Data 


Knowledge  of  wave  climatic  conditions  is  generally  required  for  ship  motion  ap¬ 
plications.  Such  applications  can  include  ship  operability  analysis  [16],  analysis  of 
ship  structural  loads  [17],  and  assessment  of  ship  capsize  probabilities  [18,  19,  20]. 

Wave  climate  data  are  most  commonly  presented  in  the  form  of  scattergrams  show¬ 
ing  joint  probability  of  occurrence  for  significant  wave  height  and  characteristic 
wave  period  (typically  peak  or  zero-crossing  wave  period).  Table  1  shows  an  ex¬ 
ample  of  a  wave  scattergram  for  the  North  Atlantic  based  on  data  from  Bales  et  al. 
[21,22]. 

Wave  climate  data  are  available  from  a  variety  of  sources.  Bales  et  al.  [21,  22] 
have  published  wave  scattergrams  developed  using  the  Spectral  Ocean  Wave  Model 
[23],  which  uses  a  hindcasting  approach  based  on  observed  wind  data.  British 
Martime  Technology  (BMT)  Global  Wave  Statistics  [24]  is  another  commonly  used 
data  source  for  wave  scattergrams,  and  is  based  on  visual  observations  from  ships. 
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Table  1:  North  Atlantic  Wave  Scattergram  Adapted  from  Bales  et  at.  [21,  22]  for 
Capsize  Risk  Analysis 
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Figure  5:  Maximum  Nominal  Wave  Steepness  Versus  Peak  Wave  Period 


McTaggart  and  de  Kat  [19]  discuss  application  of  wave  climate  data  from  these  two 
different  sources  to  prediction  of  ship  capsize  probabilities  in  the  North  Atlantic, 
and  suggest  that  the  BMT  data  give  unrealistically  large  wave  steepnesses.  This 
observation  is  based  on  comparison  with  results  from  Buckley  [25],  who  used  wave 
buoy  observations  to  develop  limiting  envelopes  of  significant  wave  height  versus 
peak  wave  period.  Buckley’s  analysis  indicates  that  nominal  wave  steepness  has  a 
limiting  value  of  0.049,  with  nominal  wave  steepness  being  given  by  the  following 
for  deep  water: 


H/X 


2  7T  Hs 

JT? 


(56) 


When  comparing  data  from  BMT  with  Buckley’s  limiting  envelope,  it  should  be 
noted  that  BMT  uses  zero-crossing  wave  period,  which  can  be  converted  to  an 
equivalent  peak  wave  period  using  Equation  (44).  Figure  5  gives  maximum  nom¬ 
inal  wave  steepness  versus  peak  wave  period  using  using  data  from  BMT,  Bales 
et  al.,  and  Buckley’s  limiting  envelope  for  the  Northern  Hemisphere.  The  nominal 
wave  steepnesses  based  on  the  BMT  data  are  much  steeper  than  Buckley’s  limiting 
envelope.  In  contrast,  the  North  Atlantic  data  from  Bales  give  steepnesses  which 
are  consistent  with  Buckley’s  envelope,  with  the  exception  of  a  single  observation 
of  Hs  =  8.5  m  and  Tp  =  9.7  s.  The  adapted  scattergram  in  Table  1  is  based  on  Bales 
et  al.,  with  the  Hs  =  8.5  m,  Tp  =  9.7  s  observation  changed  to  Hs  =  6.5  m,  Tp  = 
9.7  s. 
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6  Improved  Treatment  of  Nonlinear 
Buoyancy  and  Incident  Wave  Forces 


Reference  2  discussed  inclusion  of  nonlinear  forces  from  buoyancy  and  incident 
waves  when  predicting  ship  motions.  These  nonlinear  forces  can  be  evaluated  using 
direct  integrations  of  pressures  over  the  instantaneous  wetted  hull  surface.  The 
forces  are  initially  evaluated  using  the  ship  referenced  coordinate  system  shown  in 
Figure  6,  and  then  transformed  to  translating  earth  axes  for  solution  of  ship  motions. 


Figure  6:  Ship  Referenced  Coordinate  System  for  Large  Angular  Motions 

Reference  2  gave  equations  for  rotation  of  translation  forces  (surge,  sway,  and 
heave)  from  ship-based  to  translating  earth  axis  systems;  however,  the  transforma¬ 
tion  of  rotational  moments  (roll,  pitch,  and  yaw)  was  neglected.  For  large  amplitude 
ship  motions  given  in  translating  earth  axes,  the  order  of  coordinate  transformation 
is  yaw,  pitch,  and  roll  (i.e.,  pitch  is  given  relative  to  the  instantaneous  ship  heading, 
and  roll  is  given  relative  to  the  instantaneous  heading  and  pitch).  The  complete 
equations  for  transforming  forces  from  ship-based  to  translating  earth  axis  systems 
are  then  given  by: 

Fi  (t)  =  (t)  cos  g 6  cos  g5  —  (t)  sin  t/6  cos  g 5  +  F%  ( t )  sin  g5  (57) 

F2(t)  =  Ff(t)  (cos 7/6  cos  7/5  sin  7/4  +  sin  g6  cos 7/4) 

+  F2s(t)  (—sin 7/6  sin 7/5  sin 7/4  +cost/6  cos 7/4) 
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—  F.S(t)  sin  774  cos  775  (58) 

F3(t)  =  Ff(t)  (—cos t76  sin 775  cos 774  +  sin 770  sin 774) 

+  F2  ( t )  (sin  776  sin  775  cos  '774  +  cos  776  sin  774) 

+  F35(t)  cos  774  cos  775  (59) 

^4(t)  =  if(t)  (60) 

i^(t)  =  F5s(t)  cos  774  -  Fq  (t)  sin  774  (61) 

F6(t)  =  F55(f)  sin  774  +  F6s(t)sin774  (62) 


where  Fj  is  the  force  for  mode  j  in  translating  earth  axes,  and  Fys  is  the  force  for 
mode  j  in  ship-based  earth  axes. 

7  Numerical  Implementation  of  Seaway 
Models 


The  seaway  models  described  in  Sections  2  and  3  have  been  implemented  using  the 
Python  programming  language  [26,  27]  and  the  Numeric  Python  library  [28].  Four 
Python  classes  have  been  developed  to  model  the  following: 

•  regular  seaway  in  fixed  axes, 

•  regular  seaway  in  translating  earth  axes, 

•  random  seaway  in  fixed  axes, 

•  random  seaway  in  translating  earth  axes. 

Each  class  includes  methods  in  both  the  time  and  frequency  domains  for  comput¬ 
ing  wave  elevation,  and  wave  potential  and  its  derivatives.  To  facilitate  efficient 
computation,  the  classes  exploit  the  universal  function  capability  of  the  Numeric 
library,  whereby  input  coordinates  can  be  provided  for  a  single  location  or  an  array 
of  locations  when  evaluating  wave  elevation  or  wave  potential. 

When  computing  ship  motions  using  nonlinear  wave  excitation  and  buoyancy  forces, 
it  is  necessary  at  each  time  step  to  determine  whether  each  panel  describing  the  ship 
hull  is  below  the  water.  In  addition,  the  pressure  on  each  wetted  panel  must  be  de¬ 
termined.  These  computations  require  large  numbers  of  evaluations  of  sine  and 
cosine  functions.  In  an  attempt  to  obtain  improved  speed  performance,  functions 
were  developed  to  evaluate  sine  and  cosine  functions  using  look-up  tables.  These 
functions  were  developed  in  C  and  placed  in  a  dynamically  linked  library  (DLL) 
which  could  be  accessed  directly  by  Python.  Unfortunately,  the  new  functions  did 
not  give  any  performance  improvement  relative  to  the  Numeric  library. 
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New  Python  functions  have  been  developed  to  create  random  unidirectional  and 
multidirectional  seaways  based  on  input  wave  spectra.  To  improve  computational 
efficiency  of  time  domain  computations,  these  functions  remove  wave  components 
with  negligible  fractions  of  wave  energy  (e.g.,  less  than  one  millionth  of  total  wave 
energy.).  When  wave  components  are  removed,  remaining  wave  component  am¬ 
plitudes  are  multiplied  by  a  correction  factor  to  maintain  the  original  energy  of  the 
seaway. 

8  Computational  Requirements  for 
Simulating  Frigate  Motions  in  Random 
Seaways 


The  computational  time  required  to  predict  ship  motions  in  a  seaway  is  very  impor¬ 
tant  from  a  practical  viewpoint.  For  training  applications,  ship  motion  predictions 
should  run  faster  than  real  time.  For  engineering  analysis,  it  is  desirable  for  ship 
motion  predictions  to  run  significantly  faster  than  real-time  to  permit  examination 
of  large  numbers  of  scenarios. 

This  section  presents  some  results  of  simulations  for  DND’s  HALIFAX  class,  with 
main  particulars  given  in  Table  2.  The  same  loading  condition  was  used  for  sim¬ 
ulations  in  Reference  2.  The  hull  panelling  meshes  are  also  the  same  as  for  Ref¬ 
erence  2.  Table  3  gives  the  number  of  panels  used  for  the  three  different  meshes, 
denoted  coarse,  medium  and  fine.  Figures  7  to  9  show  all  panels  used  to  represent 
the  ship,  including  portions  below  and  above  the  calm  waterline.  The  colour  of 
each  panel  is  a  function  of  elevation  relative  to  the  calm  waterline. 

To  ensure  coursekeeping  during  simulations,  a  supplementary  yaw  stiffness  of  2  x 
109  Nm/rad  is  used,  as  was  done  in  Reference  2.  Note  the  present  work  does  not 
consider  forces  from  appendages  such  as  a  skeg  and  rudder,  which  will  contribute 
to  yaw  stiffness.  Supplementary  roll  damping  of  2  x  107  Nm/(rad/s)  represents  roll 
damping  from  viscous  forces  and  appendage  lift,  as  was  done  in  Reference  2. 

Initial  simulations  for  HALIFAX  in  random  seaways  revealed  excessive  surge  and 
sway  motions  when  using  nonlinear  terms  for  incident  wave  forces  and  buoyancy 
forces.  These  excessive  motions  were  due  to  nonlinear  wave  drift  forces,  the  lack  of 
stiffness  terms  for  surge  and  sway,  and  wave  radiation  damping  approaching  zero  as 
wave  encounter  frequency  approaches  zero.  Wave  drift  forces  arise  from  non-zero 
mean  force  terms  that  occur  when  nonlinearities  are  considered,  and  have  been  the 
subject  of  extensive  research  [29,  30,  31,  32].  Chakrabarti  [4]  provides  a  useful 
overview  of  wave  drift  forces.  To  suppress  unrealistic  surge  and  sway  motions  aris- 
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Table  2:  Main  Particulars  for  HALIFAX  Class  Frigate,  CPF  Hydroelastic  Model 
Deep  Departure  Condition 


Length,  L 
Beam,  B 

Midships  draft,  Tmid 
Trim  by  stern,  ts 
Displacement,  A 
Vertical  centre  of  gravity,  KG 
Dry  roll  radius  of  gyration  rxx 
Dry  pitch  radius  of  gyration  ryy 
Dry  yaw  radius  of  gyration  rzz 


124.5  m 

14.8  m 
4.97  m 
-0.04  m 

4655  tonnes  (fresh  water) 
6.26  m 
5.82  m 

28.8  m 
28.8  m 


Table  3:  Numbers  of  Panels  for  Different  Panel  Meshes 


Coarse 

Medium 

Fine 

Nominal  panel  size  (m2) 

5.0 

2.5 

1.0 

Number  of  panels  on  wet  hull 

448 

866 

2080 

Number  of  panels  on  dry  hull 

824 

1646 

3882 

Total  number  of  panels 

1272 

2512 

5962 
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Figure  7:  Coarse  Mesh  for  Hull  of  CPF  Hydroelastic  Model 


Figure  8:  Medium  Mesh  for  Hull  of  CPF  Hydroelastic  Model 
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Figure  9:  Fine  Mesh  for  Hull  of  CPF  Hydroelastic  Model 

ing  from  nonlinear  wave  drift  forces,  supplementary  stiffness  and  damping  terms 
were  used.  Finn  et  al.  [33]  found  it  necessary  to  use  this  approach  in  their  recent 
work  on  ship  motions  in  seaways.  When  determining  suitable  values  for  supple¬ 
mentary  stiffness  and  damping,  it  was  postulated  that  a  natural  period  of  approxi¬ 
mately  20  s  and  a  damping  ratio  (relative  to  critical  damping)  of  approximately  0.3 
would  be  appropriate  for  surge  and  sway  modes.  For  surge,  the  selected  stiffness 
was  5  x  105  N/m  and  the  selected  damping  coefficient  was  1  x  106  N/(m/s).  For 
sway,  the  selected  stiffness  was  1  x  106  N/m  and  the  selected  damping  coefficient 
was  2  x  106  N/(m/s).  The  differences  between  the  values  for  surge  and  sway  arise 
from  surge  having  negligible  added  mass  and  sway  having  added  mass  approaching 
the  mass  of  the  ship.  These  selected  supplementary  stiffness  and  damping  values 
gave  simulation  results  which  appeared  to  be  credible.  For  an  actual  ship,  stiffness 
and  damping  for  surge  and  sway  could  be  provided  by  control  systems  influencing 
propeller  speed  and  rudder  angle  to  regulate  ship  speed  and  course. 

8.1  Maximum  Recommended  Time  Step  for 
Accurate  Simulations  with  HALIFAX 

The  CPU  time  required  for  simulating  ship  motions  is  inversely  proportional  to  the 
time  step  used  during  simulations.  Consequently,  it  is  important  to  determine  the 
maximum  time  step  that  can  used  while  still  maintaining  accuracy  in  simulations. 
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Simulations  have  been  conducted  to  determine  the  maximum  time  step  that  should 
be  used  for  the  frigate  HALIFAX.  To  ensure  that  the  determined  maximum  rec¬ 
ommended  time  step  is  conservative,  the  simulations  have  been  performed  under 
conditions  whereby  results  will  be  most  sensitive  to  time  step  size.  Both  buoyancy 
and  incident  wave  excitation  forces  are  evaluated  using  the  instantaneous  water¬ 
line  at  each  time  step.  Wheeler  stretching  (Section  2.6)  is  used  for  evaluation  of 
wave  kinematic  terms.  Random  seaways  are  modelled  using  Bretschneider  spec¬ 
tra  for  significant  wave  heights  ranging  from  1  to  10  m.  The  peak  wave  period 
of  each  seaway  is  based  on  the  limiting  nominal  wave  steepness  of  0.049  deter¬ 
mined  by  Buckley  [25].  Each  modelled  seaway  has  19  wave  components,  which 
have  frequencies  of  0.2,  0.3,  . . .,  2.0  rad/s.  HALIFAX  is  travelling  at  a  moderately 
high  speed  of  20  knots  in  bow  and  stem  quartering  seas,  and  is  modelled  using  the 
medium  mesh  of  Figure  8. 

The  predicted  RMS  motions  in  Figures  10  and  11  show  that  time  step  sizes  of  0.1 
and  0.2  s  give  practically  identical  results.  A  time  step  size  of  0.5  s  gives  noticeable 
differences  in  bow  quartering  seas,  in  which  encounter  frequencies  are  greater  than 
for  stern  quartering  seas.  It  is  recommended  that  a  maximum  time  step  size  of  0.2  s 
be  used  for  simulation  of  HALIFAX  motions  in  waves. 
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-  At  =  0.1  s 

-  At  =  0.2  s 


Peak  Wave  Period  Tp  (s) 


.  At  =  0.5  s 


Peak  Wave  Period  Tp  (s) 


Figure  10:  RMS  Motions  for  HALIFAX  from  Nonlinear  Time  Domain  Computations 
with  Different  Time  Steps,  Limiting  Wave  Steepness  H/\=  0.049,  20  knots,  Stern 
Quartering  Seas  at  30  degrees 
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At  =  0.1  s 
At  =  0.2  s 


At  =  0.5  s 


Peak  Wave  Period  Tp  (s)  Peak  Wave  Period  Tp  (s) 


Significant  Wave  Height  Hs  (m)  Significant  Wave  Height  Hs  (m) 


Significant  Wave  Height  Hs  (m)  Significant  Wave  Height  Hs  (m) 

Figure  11:  RMS  Motions  for  HALIFAX  from  Nonlinear  Time  Domain  Computations 
with  Different  Time  Steps,  Limiting  Wave  Steepness  H/\=  0.049,  20  knots,  Bow 
Quartering  Seas  at  150  degrees 


24 


DRDC  Atlantic  TM  2003-190 


8.2  Influence  of  Mesh  Size  on  RMS  Motions  for 
HALIFAX 

When  computing  ship  motions  using  nonlinear  forces  from  incident  waves  and 
buoyancy,  the  computation  time  is  approximately  proportional  to  the  number  of 
panels  used  to  represent  the  ship  hull.  Consequently,  it  is  important  to  know  how 
many  panels  are  required  for  accurate  computation  of  ship  motions.  Figures  12  and 
13  show  RMS  motions  for  HALIFAX  in  steep  waves  based  on  predictions  using  the 
three  different  meshes  from  Figures  7,  8,  and  9.  The  predictions  are  based  on  a  time 
step  size  of  0.2  s.  Predictions  using  the  medium  and  fine  meshes  give  effectively 
identical  results.  Motions  predicted  using  the  coarse  mesh  give  slightly  different 
results.  The  medium  mesh  appears  to  give  the  most  suitable  balance  between  accu¬ 
racy  and  computation  time. 
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Fine  mesh 
Medium  mesh 


Peak  Wave  Period  Tp  (s) 


Coarse  mesh 


Peak  Wave  Period  Tp  (s) 


Figure  12:  RMS  Motions  for  HALIFAX  from  Nonlinear  Time  Domain  Computations 
with  Different  Meshes,  Limiting  Wave  Steepness  H/\=  0.049 ,  20  knots,  Stern 
Quartering  Seas  at  30  degrees 


26 


DRDC  Atlantic  TM  2003-190 


Fine  mesh 
Medium  mesh 


Coarse  mesh 


Peak  Wave  Period  Tp  (s)  Peak  Wave  Period  Tp  (s) 


Significant  Wave  Height  Hs  (m)  Significant  Wave  Height  Hs  (m) 


Significant  Wave  Height  Hs  (m)  Significant  Wave  Height  Hs  (m) 

Figure  13:  RMS  Motions  for  HALIFAX  from  Nonlinear  Time  Domain  Computations 
with  Different  Meshes,  Limiting  Wave  Steepness  H/\=  0.049 ,  20  knots,  Bow 
Quartering  Seas  at  150  degrees 
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8.3  Influence  of  Nonlinear  Forces  from  Buoyancy 
and  Incident  Waves  on  HALIFAX  RMS  Motions 

Computations  using  nonlinear  forces  from  buoyancy  and  incident  waves  require 
significantly  more  time  than  computations  using  linear  forces;  thus,  it  is  important 
to  know  under  what  conditions  linear  forces  can  be  used.  Linear  and  nonlinear 
computations  were  performed  using  the  HALIFAX  medium  mesh  and  a  time  step 
of  0.2  s. 

Figures  14  and  15  show  RMS  motions  for  HALIFAX  at  20  knots  in  stern  and  bow 
quartering  seas  in  moderate  wave  steepnesses  {H / A=  0.021).  Force  nonlinearities 
have  little  influence  on  ship  motions  in  moderate  wave  steepnesses,  with  the  ex¬ 
ception  of  roll  in  bow  quartering  seas  at  significant  wave  heights  greater  than  3  m. 
Figures  16  and  17  show  the  influence  of  nonlinear  buoyancy  and  incident  forces  at 
the  limiting  wave  steepness  H/X=  0.049.  Even  in  steep  waves,  force  nonlinearities 
have  little  influence  on  ship  motions,  with  the  exceptions  of  roll  and  yaw  in  bow 
quartering  seas. 
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-  Nonlinear  buoyancy  and  incident  wave  forces 

-  Linear  buoyancy  and  incident  wave  forces 

Peak  Wave  Period  Tp  (s)  Peak  Wave  Period  Tp  (s) 


Figure  14:  RMS  Motions  for  HALIFAX  from  Linear  and  Nonlinear  Force 
Computations,  Typical  Wave  Steepness  H/X=  0.021,  20  knots,  Stern  Quartering 
Seas  at  30  degrees 
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-  Nonlinear  buoyancy  and  incident  wave  forces 

-  Linear  buoyancy  and  incident  wave  forces 

Peak  Wave  Period  Tp  (s)  Peak  Wave  Period  Tp  (s) 


Figure  15:  RMS  Motions  for  HALIFAX  from  Linear  and  Nonlinear  Force 
Computations,  Typical  Wave  Steepness  H/X=  0.021,  20  knots,  Bow  Quartering 
Seas  at  150  degrees 
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-  Nonlinear  buoyancy  and  incident  wave  forces 

-  Linear  buoyancy  and  incident  wave  forces 

Peak  Wave  Period  Tp  (s)  Peak  Wave  Period  Tp  (s) 


Figure  16:  RMS  Motions  for  HALIFAX  from  Linear  and  Nonlinear  Force 
Computations,  Limiting  Wave  Steepness  H/X=  0.049,  20  knots,  Stern  Quartering 
Seas  at  30  degrees 
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-  Nonlinear  buoyancy  and  incident  wave  forces 

-  Linear  buoyancy  and  incident  wave  forces 

Peak  Wave  Period  Tp  (s)  Peak  Wave  Period  Tp  (s) 


Figure  17:  RMS  Motions  for  HALIFAX  from  Linear  and  Nonlinear  Force 
Computations,  Limiting  Wave  Steepness  H/X=  0.049,  20  knots,  Bow  Quartering 
Seas  at  150  degrees 
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8.4  CPU  Requirements  for  HALIFAX  in  Random 
Seaways 

Required  CPU  time  has  been  examined  for  HALIFAX  operating  in  random  sea¬ 
ways.  Simulations  were  conducted  in  NATO  sea  state  7  ( Hs  =  7.5  m,  Tp  =  15.0  s)  in 
both  unidirectional  waves  and  directional  waves  modelled  using  a  cosine  squared 
spreading  function  with  a  spreading  angle  of  90  degrees.  For  the  unidirectional 
waves,  the  seaway  was  modelled  using  19  wave  components  (0.2,  0.3, . . .,  2.0  rad/s) 
and  37  wave  components  (0.2,  0.25,  . . .,  2.0  rad/s).  For  multidirectional  waves,  19 
wave  frequencies  and  9  wave  directions  relative  to  the  mean  wave  direction  (-80, 
-60,  . . .,  80  degrees)  were  used.  In  the  multidirectional  seaway,  combinations  of 
wave  frequency  and  sea  direction  with  less  than  one  millionth  of  the  total  wave 
energy  were  considered  negligible  and  discarded,  leaving  a  total  of  160  wave  com¬ 
ponents.  A  time  step  size  of  0.2  s  was  used  for  the  simulations. 

Table  4  shows  the  ratio  of  CPU  time  to  simulated  time  for  the  various  cases.  The 
simulations  using  linear  buoyancy  and  incident  wave  forces  run  50  times  faster  than 
real  time  in  unidirectional  waves.  In  multidirectional  waves,  the  linear  simulations 
run  12  times  faster  than  real  time.  The  nonlinear  simulations  using  19  wave  com¬ 
ponents  for  a  unidirectional  seaway  run  4  times  faster  than  real  time. 

Table  4:  Ratio  of  CPU  Time  to  Simulated  Time  for  HALIFAX  at  20  knots,  Bow 
Quartering  Seas  at  150  degrees,  Sea  State  7,  Hs  =  7.5  m,  Tp  =  15.0  s 


Linear  or  nonlinear  buoyancy 

and  incident  wave  forces 

Number  of  wave 

components 

CPU  time/simulated  time 

Unidirectional  waves 

Linear 

19 

0.019 

Nonlinear 

19 

0.254 

Linear 

37 

0.027 

Nonlinear 

37 

0.435 

Multidirectional  waves 

Linear 

160 

0.084 

Nonlinear 

160 

1.636 

It  is  useful  to  consider  the  above  results  in  the  context  of  various  simulation  appli¬ 
cations.  Most  ship  operations,  such  as  replenishment  at  sea  and  helicopter  land¬ 
ing,  take  place  in  significant  wave  heights  less  than  7  m.  Under  such  conditions, 
differences  between  linear  and  nonlinear  ship  motion  predictions  are  very  small; 
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thus,  linear  motion  simulations  can  be  used  with  negligible  influence  on  accuracy. 
In  higher  sea  states  (i.e.,  significant  wave  heights  of  7  m  and  greater),  nonlinear 
effects  can  be  significant  and  should  be  included  in  simulations  for  which  high  fi¬ 
delity  is  required.  Examples  of  applications  requiring  nonlinear  predictions  include 
ship  capsize,  sea  loads  in  high  sea  states,  and  helicopter  securing  in  high  sea  states. 
Fortunately,  directional  scattering  decreases  in  higher  sea  states,  and  unidirectional 
seaways  can  typically  be  assumed  when  nonlinear  motion  predictions  are  required. 
It  is  likely  that  19  wave  components  are  sufficient  for  most  ship  applications  in  uni¬ 
directional  seaways.  Given  the  results  in  Table  4,  ship  motion  predictions  for  most 
applications  will  run  at  least  4  times  faster  than  real  time. 

9  Influence  of  Regular  Wave  Modelling  on 
Ship  Motions 


To  illustrate  the  influence  of  regular  wave  modelling  on  ship  motions,  numerical 
predictions  using  sinusoidal  linear  waves  and  Stokes  second  order  waves  have  been 
compared  with  experimental  results  for  the  CPF  hydroelastic  model  [34]  in  head 
seas.  Frequency  and  time  domain  predictions  use  the  medium  mesh  (Figure  8). 
The  time  domain  predictions  for  both  sinusoidal  waves  and  Stokes  second  order 
waves  include  nonlinear  buoyancy  and  incident  wave  forces.  Wheeler  stretching 
(Section  2.6)  is  applied  when  predicting  nonlinear  incident  wave  forces.  Figures  18 
and  19  show  that  the  inclusion  of  second  order  wave  effects  has  negligible  influence 
on  the  motions  of  the  CPF  hydroelastic  model  in  head  seas. 

Figures  18  and  19  show  very  close  agreement  between  linear  frequency  domain 
predictions  and  time  domain  predictions  including  nonlinear  forces  from  buoyancy 
and  incident  waves.  This  result  differs  from  previous  comparisons  given  in  Ref¬ 
erence  2.  The  differences  between  nonlinear  time  domain  predictions  from  the 
present  and  previous  studies  are  due  to  the  introduction  of  Wheeler  stretching  in 
the  present  work. 
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Experiments 
•  H/X  =  1/30 
A  H/X  =  1/20 
□  H/X  =  1/15 


Predictions,  H/X  =  1/20 
Time  domain,  second  order  Stokes  waves 
Time  domain,  sinusoidal  waves 
Frequency  domain 
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Figure  18:  Heave  and  Pitch  for  CPF  Hydroelastic  Model  in  Head  Seas,  Fn  =  0. 12 


DRDC  Atlantic  TM  2003-1 90 


35 


Experiments 
•  H/X  =  1/30 
A  H/X  =  1/20 
□  H/X  =  1/15 


Predictions,  H/X  =  1/20 
Time  domain,  second  order  Stokes  waves 
Time  domain,  sinusoidal  waves 
Frequency  domain 


Non-dimensional  Wave  Frequency:  ajj^/L/g 


Pitch 

Amplitude 

i]5/ka 


Figure  19:  Heave  and  Pitch  for  CPF  Hydroelastic  Model  in  Head  Seas,  Fn  =  0.20 
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10  Conclusions 


Models  have  been  developed  for  simulating  regular  and  random  seaways  in  both 
fixed  and  translating  axis  systems.  Translating  axis  systems  are  commonly  used  for 
evaluating  wave  forces  on  ships.  The  presented  phase  relationship  between  fixed 
and  translating  axis  systems  can  be  used  to  ensure  consistency  of  computations  and 
graphically  presented  results.  For  regular  waves,  the  nonlinear  wave  surface  can 
be  modelled  using  Stokes  second  order  theory.  The  influence  of  nonlinearities  on 
wave  kinematics  can  be  modelled  using  Wheeler  stretching,  which  can  be  applied 
to  both  regular  and  random  seaways. 

Commonly  used  wave  spectral  models  have  been  presented.  Directional  wave  ef¬ 
fects  can  be  modelled  by  applying  a  cosine  squared  spreading  function  to  a  point 
wave  spectrum,  or  by  a  ten  parameter  spectrum,  which  includes  directional  terms. 

A  brief  review  of  wave  climate  data  has  discussed  available  sources.  When  select¬ 
ing  significant  wave  height  and  peak  wave  period  for  modelling  random  seaways, 
care  should  be  taken  to  not  exceed  the  limiting  nominal  wave  steepness  of  0.049 
based  on  wave  buoy  observations. 

A  parametric  investigation  was  performed  to  investigate  sensitivity  of  RMS  mo¬ 
tions  to  the  following  parameters  that  significantly  influence  required  computational 
time:  hull  panel  size,  time  step  size,  and  usage  of  linear  or  nonlinear  forces  from 
buoyancy  and  incident  waves.  Motions  were  predicted  for  HALIFAX  at  20  knots  in 
stern  quartering  and  bow  quartering  random  seaways.  Computations  using  different 
panels  sizes  show  convergence  of  RMS  motions  when  using  a  medium  mesh  with  a 
nominal  panel  area  of  2.5  m2,  giving  1646  panels  on  the  wet  portion  of  the  hull  and 
2512  panels  on  the  dry  portion  of  the  hull.  Computations  using  different  time  step 
sizes  show  good  convergence  of  RMS  motions  for  a  time  step  of  0.2  s.  Inclusion 
of  nonlinear  forces  due  to  buoyancy  and  incident  waves  significantly  influences  roll 
and  yaw  predictions  in  bow  quartering  seas.  Nonlinear  effects  are  small  for  other 
modes  in  bow  quartering  seas  and  for  all  modes  in  stem  quartering  seas. 

For  simulations  in  random  seaways,  the  ratio  of  CPU  time  to  simulated  time  is 
highly  dependent  upon  the  number  of  wave  components  and  whether  nonlinear 
forces  from  buoyancy  and  incident  waves  are  computed  at  each  time  step.  Using 
a  800  MHz  Pentium  III  computer,  the  ratio  of  CPU  time  to  simulated  time  ranged 
from  0.02  to  1.6. 

Comparisons  of  heave  and  pitch  predictions  for  HALIFAX  with  experimental  re¬ 
sults  in  regular  waves  show  very  good  agreement,  with  nonlinear  effects  being 
small. 
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spectral  normalization  factor 
wave  amplitude 

amplitude  of  wave  component  i 
beam 

British  Maritime  Technology 
dynamically  linked  library 
Froude  number 
directional  spreading  function 
gravitational  acceleration 
nominal  wave  steepness 
significant  wave  height 

significant  wave  height  of  bimodal  spectrum  component  i 

height  of  centre  of  gravity  above  baseline 

incident  wavenumber 

incident  wavenumber  of  wave  component  i 

ship  length  between  perpendiculars 

directional  spreading  function  of  bimodal  spectral  component  i 

V  th  moment  of  spectrum 

number  of  incident  wave  components 

spectral  directional  spreading  parameter 

root-mean- square 

roll  radius  of  gyration 

pitch  radius  of  gyration 

yaw  radius  of  gyration 

wave  spectral  density 

draft  at  midships 

peak  wave  period 

zero-crossing  period  of  spectrum 

average  wave  period 

trim  by  stern 

reference  time 

ship  forward  speed 

horizontal  plane  coordinates  in  translating  earth  axes 
horizontal  plane  coordinates  in  fixed  axes 
vertical  coordinate  relative  to  calm  water  surface 
normalization  factor  for  JONS  WAP  spectrum 
wave  direction  relative  to  ship 
peak  enhancement  factor  for  JONSWAP  spectrum 
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wave  spreading  angle 

frequency  increment  of  wave  component  i 

phase  lead  of  incident  wave  in  translating  earth  coordinate  system 

phase  lead  of  incident  wave  in  fixed  coordinate  system 

phase  lead  of  incident  wave  component  i 

incident  wave  elevation 

motion  displacement  in  mode  j 

exponent  for  JONS  WAP  spectrum 

spectral  shape  parameter  of  bimodal  spectrum  component  i 

wave  heading  (from)  in  fixed  coordinate  system 

mean  wave  heading  of  seaway  (from)  in  fixed  coordinate  system 

water  density 

standard  deviation 

incident  wave  potential  in  time  domain 

incident  wave  potential  in  frequency  domain 

ship  heading  (toward)  in  fixed  coordinate  system 

incident  wave  frequency 

incident  frequency  of  wave  component  % 

wave  encounter  frequency 

peak  wave  frequency 

peak  wave  frequency  of  bimodal  spectrum  component  i 
ship  mass  displacement 
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of  the  location  of  the  free  surface  on  wave  kinematics  can  be  modelled  using 
Wheeler  stretching.  Commonly  used  wave  spectra  and  sources  of  wave  climate 
data  are  presented  as  reference  information.  Sample  computations  for  the  frigate 
FIALIFAX  illustrate  the  influence  of  simulation  parameters  on  predicted  RMS  mo¬ 
tions  and  required  computational  time.  Required  computational  time  increases 
significantly  when  including  nonlinear  forces  from  buoyancy  and  incident  waves. 
In  most  cases,  simulations  run  faster  than  real  time. 


14.  KEYWORDS,  DESCRIPTORS  or  IDENTIFIERS  (technically  meaningful  terms  or  short  phrases  that  characterize 
a  document  and  could  be  helpful  in  cataloguing  the  document.  They  should  be  selected  so  that  no  security  classification 
is  required.  Identifiers,  such  as  equipment  model  designation,  trade  name,  military  project  code  name,  geographic 
location  may  also  be  included.  If  possible  keywords  should  be  selected  from  a  published  thesaurus,  e.g.  Thesaurus  of 
Engineering  and  Scientific  Terms  (TEST)  and  that  thesaurus-identified.  If  it  not  possible  to  select  indexing  terms  which 
are  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title). 
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